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ABSTRACT 

We quantify the effects of electron thermal conduction on the properties of hot accretion flows, 
under the assumption of spherical symmetry. Electron heat conduction is important for low accretion 
rate systems where the electron cooling time is longer than the conduction time of the plasma, such 
as Sgr A* in the Galactic Center. For accretion flows with density profiles similar to the Bondi 
solution (n[r] oc r~ 3 / 2 ), we show that heat conduction leads to super- virial temperatures, implying 
that conduction significantly modifies the dynamics of the accretion flow. We then self-consistently 
solve for the dynamics of spherical accretion in the presence of saturated conduction and electron 
heating. We find that the accretion rate onto the central object can be reduced by ~ 1 — 3 orders 
of magnitude relative to the canonical Bondi rate. Electron conduction may thus be an important 
ingredient in explaining the low radiative efficiencies and low accretion rates inferred from observations 
of low- luminosity galactic nuclei. The solutions presented in this paper may also describe the nonlinear 
saturation of the magnetothermal instability in hot accretion flows. 
Subject headings: accretion, accretion disks - Galaxy: center 



1. INTRODUCTION 

The conversion of gravitational binding energy of ac- 
creting matter into radiation is thought to account for the 
high lumino sities observed fro m many compact objects 
(e.g., AGN: lLvnden-Bell|[l969l ). Most of the time, how- 
ever, systems that appear to harbor compact objects (low 
luminosity AGN, X-ray binary systems in quiescence) do 
not radiate at the levels one would expect from simple 
estimates of the efficiency at which binding energy is con- 
verted to radiation. This implies that these systems have 
low radiative efficiencies and/or low accretion rates, and 
that the accreting gas may thus retain most of its binding 
energy in the form of heat. The ultimate energy sink for 
such hot radiatively inefficient accretion flows (RIAFs) 
remains an area of active research. T he energy may 
be pr imarily advected in with the flow (|Naravan fc Yil 
1994), or it may be directed into driving bulk convec- 
tive m otion (|Quataert fe Gruzinovll2000al: iNaravan et al.1 
I2000D or a global outflow (jBlandford fc Begelmanlll999f) . 
Numerical simulations favor th e latte r possibilities (e.g., 
Igumenshchev fc Abramowiczl Il999t IStone fc Pringld 
20011: lHawlev fc Balbudl2002T . 



Another possible outlet for the accretion energy is a 
heat flux due to thermal conduction. In order to main- 
tain a hot accretion flow, the timescale for electrons and 
protons to exchange energy by Coulomb collisions ex- 
ceeds the timescale for the plasma to flow into the central 
black hole (e.g., Rees et al. 1982). Thus the inflowing 
plasma is collisionless and energy transport by heat con- 
duction may be dynamically important. Direct observa- 
tional estimates of the electron mean-free path on scales 
resolved by Chandra in nearby galactic nuclei - includ- 
ing our Galactic Center - give values ~ 0.02 — 1.3 times 
the gravitational radius of influence of the central black 
hole (~ 10 5 — 10 6 Schwarzschild radii). Thus the plasma 
is relativel y collisionless even at large distances fro m the 
black hole (|QuataertJl200l iTanaka fe Menoull2006f ) . Un- 
der these conditions, the magnitude of the heat flux is 
uncertain and depends on the geometry of the magnetic 



field and the effective collision rate due to plasma waves 
and instabilities. It is plausible, however, that thermal 
conduction will occur at close to its maximal saturated 
rate, with a heat flux F oc pv^, w here p is the density an d 
v e the electron thermal speed (jCowie fcMcKedll977h 
This appears to be true in the solar wind (| Salem et al.1 
2003) , a well-studied example of a macroscopically colli- 
sionless plasma. For non-relativistic electrons, the ratio 
of the inflow time to the electron conduction time in a 
hot accretion flow may be as large as ~ v e /v r (see equa- 
tion [Tj below), where v r is the radial flow velocity. For 
subsonic inflow (with respect to the ions), this ratio is 
> ^Jm p /m e ~ 40, highlighting the possible imp ortance 
of con duction in hot accretion flows. In addition. iBalbud 
(2000) has shown that low collisionality conducting plas- 
mas are unstable to convective-like motions in the pres- 
ence of an outwardly decreasing temperature gradient 
(the MTI or magnetothermal instability). Although the 
interplay between the MTI and the magnetorotational 
instability (MRI) is uncertain, it is plausible that the 
MTI will lead to a dynamically important heat flux in 
hot accretion flows. 

Any significant heat flux from the inner (hotter) to 
the outer (colder) parts of an accretion flow can act to 
reduce the inflow rate. T his is true if energy is trans- 
ported by ra diation fe.g.. lOstriker et alj|1976j ). viscous 
stress es (e.g.. [Blandford fc Begelmanlll999D . convection 
(e.g., iQuataert fc Gruzinovll2000a[ K or thermal conduc- 
tion, the focus of this paper. Gas in the outer region 
of a spherical accretion flow near the sonic point has a 
thermal energy comparable to its gravitational binding 
energy. Additional heating will drive the temperature of 
the gas above the virial temperature. If the gas is unable 
to cool efficiently, this increase in temperature will lead 
to a decrease in accretion rate by either decreasing the 
radius of the sonic point or by driving a thermal outflow. 

Motivated by these considerations, this paper as- 
sesses the importance of thermal conduction using a 
simple one-dimensional (spherical) model of hot accre- 
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tion fl ows with saturated conduction. iTanaka fc Menoul 
(2006) have carried out a related analysis and find 
that the accretion flow can spontaneously produce ther- 
mal outflows driven in part by conduction (as in early 
models of the solar wind). Their analysis is two- 
dimensional but self-similar in radius, whereas ours is 
one-dimensional but does not assume self-similarity; we 
also treat saturated heat conduction somewhat differ- 
ently than they do (see 52]). Tanaka & Menou's as- 
sumption of self-similarity enforces a density profile 
that varies as r -3 / 2 , whereas simulations of RIAFs 
consis tently find densit y profiles shallower than this 
(e.g., Stone ct al. 1999, Igumenshchev & Abramowicz 
19991 iStone fc Pringld 120011 lHa wlev fc Balbusl 120021 
Igumenshchev et al]|2003D . 



We begin in [}2]by outlining the energy equation in the 
presence of saturated heat conduction. For completeness 
we consider both the ion and electron energy equations, 
though our focus will be on the electrons throughout this 
paper; we also quantify when electron heat conduction is 
likely to be important as an energy transport mechanism 
relative to radiative cooling. We then present order of 
magnitude arguments that demonstrate the importance 
of heat conduction in hot accretion flows. In particular, 
we solve for the electron temperature in the presence of 
heat conduction under the simplifying assumption that 
the density profile of the accretion flow is a fixed power- 
law in radius ($3}. In $4] we relax the power-law density 
profile assumption of SJJ] an d self-consistently solve for 
the density, temperature, and radial velocity of the flow 
neglec ting angular momentum; this is a version of the 
classic Bondi (1952) accretion problem. We summarize 
and discuss the implications of our findings in ij5] 

2. ENERGY EQUATIONS 

As noted in the Introduction, in a hot accretion flow 
the timescale for electrons and ions to exchange energy by 
Coulomb collisions exceeds the timescale for the plasma 
to flow into the central black hole. Thus the flow likely 
develops a two-temperature structure, although the ex- 
tent to which the particle species are coupled through 
other processes such as wave-particle interactions is un- 
clear. In general, then, the energy of the electrons and 
ions must be tracked separately. In spherical symmetry, 
the electron and ion energy equations can be written as 

v r nT e ^- = g+ - -^--^ (r 2 Q e ) - q~ 7 (1) 



dr 



dr 



and 



v r nTi^p- = qf - (r 2 Qi 
dr r z dr 



(2) 



where v r = —dr/dt is the radial velocity (defined here 
to be positive inward), n is the particle number density 
(n e = rii = n by charge neutrality, where the subscripts 
e and i denote properties of the electron and ion fluids, 
respectively), s is the specific entropy, q + is the heating 
rate per unit volume, and q~ is the electron cooling rate 
per unit volume. 1 In equations {1} and ^ we have ne- 
glected energy exchange between ions and electrons (e.g., 
via Coulomb collisions) and ion radiative cooling. The 

1 We work in units in which Boltzmann's constant kg = 1 so that 
temperature has units of energy, although we express our numerical 
results in Kelvin. 



last term on the right-hand side of equations ([1} and |2]) 
is the divergence of the radial conductive heat flux 



Q, 



dT e 



dr 



(3) 



Because the accreting plasma is collisionless, we model 
the thermal conductivity n using a saturated thermal 
conductivity (|Cowie fc McKedli97l : 

K e ee a c v c nr, (4) 

where a c is a dimensionless parameter that measures the 
strength of the conduction and 



for v e 
for v e 



< c 
> c 



(5) 



is a function that interpolates smoothly between the elec- 
tron thermal speed at low temperatures and the speed of 
light c as the electrons become relativistic. For the ions, 
which are always non-relativistic, ki = a c Vinr where Vi 
is the ion thermal velocity. W e use a slightly differen t 
form for the saturated flux than lCowie fc McKeel (|1977f) . 
Their flux is proportional to density and temperature, 
while we take the flux to be proportional to the tem- 
perature gradient. For our problem, the temperature 
and density are b ound ary conditions at large radii; the 
iCowie fc McKeel ()1977t ) prescription would lead to the 
heat flux being specified by these boundary conditions, 
which is unphysical since we want to be able to solve 
for the heat flux rather than have it predetermined. We 
note, however, that our assumption of Q oc dT e /dr, while 
physically reasonable, cannot be rigorously derived in the 
collisionless limit, where the heat flux depends on a non- 
local integral of the electron temperature along magnetic 
field lines (see, e.g., §VII of Snyder et al. 1997). 2 

Assuming that a fraction f e (fi) of the available ac- 
cretion energy goes into heating the electrons (ions), the 
electron and ion heating rates per unit volume are given 

by 

+ f GMM 

where fi — 1 — f e . Integrating the total (electron + ion) 
volumetric heating rate over all radii yields a net heat- 
ing rate of GMM/Ri n , where Ri n is the inner radius of 
the flow. We set Ri n — 3R,s — 6GM/c 2 throughout this 
paper so that the total heating rate is « 0.16Mc 2 . Note 
that even for spherical accretion (i.e., in the absence of 
differential rotation and the MRI), a significant fraction 
of the gravitational potential energy of the inflowing gas 
may be converted into heat. In particular, even initially 
weak magnetic fields are compressed by flux freezing until 
they become strong, leading to reconnection a nd signifi- 
cant heating (Ig umenshchev fc Naravanll2002f h 

Because v e Vi unless T e < 10~ 3 X"i, electron heat 
conduction tends to be more important than ion heat 
conduction. Our focus in this paper will thus be on elec- 
tron heat conduction. We note, however, that electron 
heat conduction is important only when electron radia- 
tive losses are relatively small. Otherwise, most of the 
energy supplied to the electrons is radiated away and 

2 If wave-particle scattering limits the mean free path of particles 
to be < r at any radius, then Q oc dT e /dr is again well-justified. 
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there is little energy available to redistribute by con- 
duction. Throughout this paper we thus set q~ = 
in equation (Q]). This is valid only when the cooling 
time of the electrons in the accretion flow is long com- 
pared to the conduction time, r/v c . At the low accretion 
rates where this assumption is appropriate (see below), 
and inside ~ f 3 Rs where the electrons are relativistic, 
synch rotron cooling is the domin ant cooling mechanism 
(e.g., iQuataert fc Naravanl fT999). Synchrotron cooling 
for thermal electrons at temperature T e can be neglected 
only for 



M 



Mi 



EDD 



\ 10/ Vo.i/ vo.i 



T, 



10 n K 



3/2 



(7) 

where Medd = £edd/(0.1c 2 ) is the Eddington accretion 
rate, the radial velocity is assumed to be a fraction a of 
the free-fall velocity, the magnetic field has been normal- 
ized to the virial temperature via f3 = nTi / (B 2 / 8ir) , the 
average Lorentz factor for a thermal distribution of elec- 
trons is given by 3T e /(m e c 2 ), and v c w c for relativistic 
electrons. Expression ((TJ) is considerably more restrictive 
than the requirement that M < c* 2 A:/edd for Coulomb 
collisions to be unimportant, which is necessar y to main- 
tain a two-temperature flow (jRees et al.l fT982). There is 
therefore a large range of accretion rates where a RIAF 
is possible but where electron heat conduction is unim- 
portant because the electrons radiate away most of their 
energy. 3 Our focus here is on very low accretion-rate 
systems where electron radiative losses are small. This 
includes Sgr A* in the Galactic Center and other very 
low-luminosity AGN and X-ray binaries in quiescence. 
Ion heat conduction may be important at higher accre- 
tion rates, but this is harder to assess because the ion 
conduction time is comparable to the inflow time. 

We note that in their tre a tment of heat conduction 
in RIAFs, iTanaka fc Menoul (|2006f ) use a single energy 
equation, as we do. Interpreted as electron conduction, 
their results- like ours- are only relevant for very low 
accretion rates. Interpreted as ion conduction, their <f>g 
should be < 0.01 (and their A < 0.1) since using the 
ion thermal speed in the conductivity introduces an ad- 
ditional factor of ^m e /m p . For these values, their solu- 
tions show that ion conduction has only a modest effect 
on the properties of RIAFs. 

The appropriate values for a c and f e are set by com- 
plicated physical processes that are not fully understood. 
The electron heating parameter must satisfy f e < 1 
on energetic grounds . For reaso ns explain e d in, e.g., 
IQuataert fc Gruzinovl (|1999h and iQuataertl (|2003h . we 
suspect that the electrons likely receive a significant frac- 
tion of the available gravitational potential energy in RI- 
AFs, so that / e ~ 0.1 — 1 is reasonable. The dimcn- 
sionless conductivity a c must also satisfy a c < 1. We 
show solutions with a c > 1, however, in order to illus- 
trate the change in dynamics as the inflowing plasma 
becomes isothermal. In the Appendix we place an upper 
limit on the ion conductivity of a c < 0.1 based on pitch 
angle scattering by instabilities generated by the nonlin- 
ear evolution of the MRI in a collisionless plasma. This 

3 Note that at high accretion rates electron conduction could 
still be important at radii > 10 3 Rs where the electrons are non- 
relativistic and synchrotron losses are small. 



argues against the importance of ion conduction. The 
corresponding constraint on the electron conductivity is 
much weaker (eq. IA3|) . 

With equations ([I])-© describing the electron ener- 
getics, we now proceed to assess the importance of heat 
conduction in RIAFs by considering a couple of model 
problems. 

3. ORDER-OF-MAGNITUDE CONSIDERATIONS 

To quantitatively motivate the importance of thermal 
conduction, we consider a simple model in which a frac- 
tion f e of the accretion energy (O.lMc 2 ) is carried out 
to large radii by a saturated conductive heat flux with 
a c ~ 1, in which case 0.1 / e Mc 2 « A-Kpr 2 v 3 (m e /m p ). For 
constant M, this implies 

\ 2 / 3 / „ \ 2 / 3 

O-i/W— -ft) ~30(/e 



1± 

T 




10 5 i?c 



(8) 

where we have assumed that the electrons are non- 
relativistic, which is true at large radii, r/Rs > 10 3 . 
Equation ([8]) shows that if even a small fraction of the 
available energy is transported outwards by conduction 
(fe ^ 0.01), the heat flux can unbind the gas at large 
radii, thus significantly changing the flow dynamics. 

The above analysis assumes that M is constant and 
thus that the density profile of the flow is given by 
p oc r~ 3 / 2 , the canonical result for spherical accre- 
tion. Models of rotating radiatively inefficient ac- 
cretion flows show, however, that the dynamics of 
the flow can be quite different from tha t of spheri- 
cal accretion. In particular, convection (I Stone et all 
119991: IQuataert fc Gruzinovl l2000at iNaravan et al.ll2000D 
and/or thermally driven winds (Blandford & Bcg elmanl 
[T^IStone fc Prinddl20lnl : lHawlev~fc Balbusll2002T ) can 
strongly suppress the rate of mass accretion onto the cen- 
tral black hole. A simple model for such a flow is one 
in which the density profile is given by n oc r~ p , so that 
M cx r 3 / 2 ~ p is steady in time but varies with radius. Nu- 
merical simulations favor p w 1/2 — 1. For p < 3/2, the 
mass accretion rate decreases at small radii, and there is 
less energy available to be conducted out to large radii. 
In this case, energy transport by conduction has less of 
an effect on the dynamics of the flow. To quantify this, 
we again assume that a fraction f e of the accretion en- 
ergy at small radii (0.lM[Rs]c 2 ) is carried out to large 
radii by a saturated conductive heat flux; allowing for 
M cx r 3 ' 2 ~ p , tn e temperature at large radii is given by 

2/3 

"'■ ' ' (9) 




r 


p-l/2\ 


Rs. 





for p > 1/2 and 




Rs 



2/3 



(10) 



for p = 1/2. For p = 1/2, equation (flTJ|l shows that the 
electron temperature required to carry the heat flux out 
is always sub-virial and thus conduction is not dynami- 
cally important. More generally, from equation (|9]), and 
assuming a fiducial outer-radius of r ~ 10 5 Rs, we find 
that the electron temperature is super- virial for p > 1. 
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To expand on the above analysis, we solve the full elec- 
tron energ y equation (eg. [1]) using a standard relaxation 
algorithm ([Press et al J 11992) for a fixed density profile 
n oc r~ p . We take the the radial velocity to be given by 



H 



V<f, 



GM 



(11) 



where H ~ r is the vertical scale-height for a hot ac- 
cretion flow, = rQ is the rotational velocity and 

£1 = W GM/r 3 is the Keplerian rotational frequency. 
There are two boundary conditions for this problem since 
we are specifying both the radial velocity and density 
profiles. A reasonable set of boundary conditions is to 
specify the temperature at the outer radius and to require 
that dT e /dr\n in = 0, i.e., that the conductive heat flux 
is zero at the inner radius (so that there is not a source 
of energy at the horizon). Our results are all insensitive 
to the precise inner boundary condition (for modest val- 
ues of dT e /dr\n in ). The parameters that determine the 
electron temperature in this model are a c /a, p, and f e . 

Figure Q] is a plot of T e (r) for different choices for the 
density power-law index p. The values for the other pa- 
rameters are a c /a = 1, f e = 0.1 and an outer tempera- 
ture T out — 2 x 10 7 K (the temperature observed in the 
Galactic Center; Baganoff et al. 2003). Figure [T] shows 
that larger values of p lead to larger electron temper- 
atures (see also eqs. [8]- [10]). This is because larger 
values of p correspond to larger accretion rates and elec- 
tron heating rates at small radii. In the presence of sig- 
nificant heat conduction some of the energy supplied to 
the electrons at small radii is conducted to large radii 
(for a c ~ a, the heat conduction time is ~ v c /v<p > 1 
times smaller than the inflow time; thus most of the en- 
ergy supplied to the electrons is conducted to large radii, 
as assumed in the order of magnitude estimates above). 
Larger values of p thus correspond to a larger heat flux at 
large radii, and a larger electron temperature is needed 

to carry this heat flux (eqs. [3]- [5] imply Q r oc T^ 2 at 
large radii where the electrons are non-relativistic). For 
sufficiently large values of p > 1 the solutions shown in 
Figure Q] are unphysical: the electron temperature be- 
comes super- virial, indicating that the electron pressure 
gradient is dynamically important and that the electrons 
are unbound from the system (see also eq. [9]). 

Using solutions with p — 1 and p = 1/2 as reference 
calculations (values motivated in part by simulations of 
RIAFs) , we explore how the temperature profile varies as 
a function of a c /a and f e . These calculations are shown 
in Figures [5] and [31 The trends in these figures can all 
be understood relatively easily. Larger values of f e cor- 
respond to larger electron heating rates and thus lead to 
higher electron temperatures. Larger values of a c /a, on 
the other hand, correspond to very efficient heat conduc- 
tion that drives the electron temperature profile towards 
isothermality. Note that in the case of p = 1/2, the elec- 
tron temperature is never super- virial (Fig. [3]), consistent 
with our order-of-magnitude estimate in equation (jlOp . 

4. SPHERICAL ACCRETION WITH CONDUCTION 

The fact that the electron temperature calculated in 
the previous section is sometimes super-virial implies 
that conduction can significantly influence the dynam- 
ics of the accretion flow; this motivates the investigation 



of this subsection, in which we solve for the density and 
velocity profiles self-consistently. We consider a simple 
model of a hot accretion flow, assuming spherical sym- 
metry, constant M, and neglecting angular momentum. 

With the above assumptions the basic dynamical equa- 
tions for the density and radial velocity are 



M = 47rr pv r 

dv r 1 dP e GM 
dr p dr ( r - R s ) 2 



= 0, 



(12) 
(13) 



where P e is the electron p ressure and we have used the 
iPaczvnsky fe Wiital (jl980l ) gravitational potential to en- 
sure that all the solutions have a sonic radius outside the 
horizon. In addition to equations (fT2"[) and lfT3]) we solve 
equation (Q]) for the electron temperature, which can also 
be written as 



dr 



7e 



1 



„ dn f P GMm n nv r 

-v r T e —+ J - j-H— £ 

dr r z 



a c a 

—o~r I r~nv c 
r z dr 



(14) 

where the entropy gradient has been rewritten in terms 
of density and temperature gradients and the internal en- 
ergy has been expressed in terms of an effective adiabatic 
index 



7eS - 



1 + 2v 2 Jc 2 

2 + 5v 2 /c 2 



5/3 for v e <C c 
4/3 for v e > c 



(15) 



that interpolates smoothly between nonrelativist ic and 
relatiy istic energies. Equations lfll^ - lfl"4]) are the iBondl 
(1952) accretion problem except that there is now ex- 
plicit heating and an additional dependent variable 
(dT e / dr) due to the conductiv e flux term in equation 
((14]) . We note that iGruzinovi (|1998fl has also consid- 
ered Bondi accretion in the presence of heat conduction, 
although his form of the conductive heat flux is quite 
different from ours and leads to different results. 

In general, the pressure in equation (|13() is given by 
the sum of the electron and ion gas pressures (radia- 
tion pressure is negligible since the plasma is optically 
thin and we do not include the effects of magnetic fields) 
and one needs to solve equation ([2|) for the ions in addi- 
tion to equation (fT4|) for the electrons. Since our goal is 
to understand the effects of conduction on the solution, 
however, it is reasonable to neglect the ion pressure in 
equation ([T3]) and focus solely on the electron energetics. 
The ion temperature in hot accretion flows is approxi- 
mately virial, Tj oc 1/r, where the constant of propor- 
tionality depends upon the assumption one makes about 
the fraction fi of accretion energy that goes into heating 
the ions. In Sj5[ we briefly discuss the effect of including 
the ion energy equation in our solutions. 

Eliminating the density as a dependent variable in fa- 
vor of the radial velocity, the momentum equation can 
be expressed as 



1 dv r 
v r dr 



2cj 



d4 

dr 



GM 



(r - Rs) 2 ' 



(16) 



where 



P e / p is the isothermal sound speed. Equa- 



tion (|16p is singular when the factor in parentheses on the 
left-hand side is zero. By analogy with the Bondi prob- 
lem, we expect a family of subsonic solutions as well as a 
unique solution that makes a transition between subsonic 
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and supersonic flow at an intermediate radius R son ic (the 
sonic point). The solution that passes through the sonic 
point is selected by setting both the left-hand and right- 
hand sides of equation (|16[) to zero at R SO nic- Specifying 
two additional boundary conditions over-constrains the 
problem, so that R SO nic (and therefore M) is uniquely 
determined as an eigenvalue. As in the previous section, 
the additional two boundary conditions are the tempera- 
ture at the outer radius T out and the zero-flux condition 
(dT e /dr = 0) at the inner radius. Note that the problem 
solved in this section is thus a 3 point boundary value 
problem, with boundary conditions at the inner radius, 
sonic point, and outer radius. We solve this problem with 
a relaxation algorithm. 

Figures IU [5] and [6] show the temperature, density and 
velocity profiles for a set of calculations with f e = 0.3, 
T out = 2 x 10 7 K and various values of a c . Our choice of 
T out corresponds to a gravitational sphere of influence of 
the black hole of GM/c 2 s \ Rout « 10 5 R S . As in Figures H 
and [31 large values of a c drive the electron temperature 
to be isothermal. The electron temperature is maximized 
(and the density at small radii minimized) for a c ~ 1. 
The differences relative to the standard Bondi problem 
are quantified more clearly in Figures [7] and [8j which 
show the sonic radius and accretion rate as a function 
of a c for several values of f e . In all of our calculations, 
the standard Bondi problem is recovered as a c — > and 
f e — > 0, though note that we are employing an adiabatic 
index that varies with radius so the a c — > solution does 
not correspond to a constant 7 as in the Bondi problem. 
In the limit of infinite conductivity, a c 3> 1, the electron 
temperature profile becomes isothermal, independent of 
f e , and the problem reduces to the standard Bondi prob- 
lem with 7 = 1. 

For intermediate values of a c ~ 1, however, the so- 
lution differs significantly from the canonical Bondi so- 
lution. In particular, Figure [7] shows that for moderate 
values of f e and a c ~ 0.1 — 1, the sonic point moves in 
from ~ 10 5 i?5 to ~ 10 - 100 Rs- The accretion rate 
is correspondingly reduced relative to the adiabatic or 
isothermal Bondi rate by ~ 1 — 3 orders of magnitude 
(Fig. [8j). Physically, this occurs because accretion at the 
Bondi rate would imply so much heat conduction to large 
radii that the electron temperature would be super- virial, 
as found in the previous section. To maintain T e < T V i r 
requires that the accretion rate decrease so that less en- 
ergy is conducted out to large radii. This occurs via a 
decrease in the radial velocity (Fig. [S]) relative to free 
fall and a flattening of the radial density profile relative 
to the r -3 / 2 Bondi solution (Fig. [5]). 

The fact that the accretion rate for nonzero f e differs 
from the adiabatic Bondi rate at low a c can be under- 
stood as follows. In the outer, subsonic region, the mo- 
mentum equation is dominated by hydrostatic balance: 

d (nT e ) ( GMm p n = Q 



dr 



+ 



For a c = 0, the energy equation is given by 
d ( T e 



T dn f e GMm p n _ Q 



dr \7e — 1/ dr 
Combining these two relations gives 



(18) 



d_ 

dr 



7e 



7e 



1" 



GMm p GMm v 

fe 



which for the form of 7 e we have chosen gives a quadratic 
relation for T e . For constant 7, T e oc 1/r and n oc r~ p , 
with 

P== (l + / e )(7-D- (20) 
The assumption of hydrostatic equilibrium breaks down 
for p > 3/2 (the free-fall value), or 



7 < Icrit 



5 (1 + f e ) 

3(l + 5/ B /3)' 



(21) 



0, (19) 



For 7 < j cr it a supersonic free-fall solution is obtain- 
able, while for larger 7 it is not. Since "fcrit < 5/3 for 
nonzero f e , the assumption of hydrostatic equilibrium 
holds for a larger range of radii than it would in the ab- 
sence of heating (until the temperature increases and 7 e 
decreases enough to violate condition [H]). There is a 
corresponding decrease in the sonic radius and accretion 
rate even for a c — > 0, as can be seen in Figures [7] and [51 
Also, since 7 e > 4/3, equation (|2ip indicates that only 
solutions with f e < 0.6 will be able to make a transition 
to supersonic flow as a c — > (in a Newtonian poten- 
tial). For larger values of f e the flow is in approximate 
hydrostatic equilibrium down to very small radii, where 
the singularity in the pseudo-Newtonian potential forces 
a sonic transition. 

5. DISCUSSION 

We have shown that electron heat conduction can have 
a significant effect on the properties of radiatively in- 
efficient accretion flows (RIAFs). We considered two 
model problems to demonstrate this, both of which as- 
sume spherical symmetry. First, for the canonical Bondi 
problem, we find that if a few percent of the accretion 
power is transferred to large radii by heat conduction, the 
heat flux can unbind the gas at large radii (eq. [5]). Gen- 
eralizing this calculation to allow for a power-law density 
profile with n(r) oc r~ p (Sj3]and Figs. DE]), we solved for 
the electron temperature profile in the presence of sat- 
urated heat conduction; p — 3/2 is appropriate for the 
spherical Bondi problem, while p ~ 1/2 — 1 is typically 
found in numerical simulations of rotating radiatively in- 
efficient accretion flows. Our calculations show that for 
density profiles with p > 1, and for plausible values of 
the electron conductivity and electron heating, the elec- 
tron temperature is super- virial at large radii (see Fig.[I} 
[3]). This indicates that the effects of conduction on the 
dynamics of the accretion flow must be self-consistently 
taken into account. 

Motivated by this initial calculation, we solved the 
problem of spherical accretion in the presence of heat 
conduction and electron heating. For appreciable elec- 
tron heating rates (f e ~ 0.1) and dimensionless conduc- 
tivities a c ~ 0.1 — 1 (sec 521 for the details of the electron 
energy equation), the accretion rate is reduced by ~ 1 — 3 
orders of magnitude relative to the canonical Bondi rate 
(21 Fig. [SJ. Physically, this is because some of the en- 
ergy dissipated as heat is conducted to large radii where 
the increase in pressure stifles the inflow of matter. This 
decrease in M relative to the Bondi accretion rate is con- 
sistent with observations of the Galactic Center, which 
indicate that the accretion rate is 2 — 3 orders of mag- 
nitude lower than the Bondi rate (iQuataert fc Gruzinovl 
2000B iBower et alj|2003t iMarrone et al.ll2006D . 



() 



Johnson & Quataert 



As discussed in $2] electron heat conduction is only im- 
portant at very low accretion rates < 1CP 5 — 10~ 6 MedDj 
when the cooling time of the electrons via synchrotron ra- 
diation is longer than the conduction time in the plasma 
(eq. [7]). Such accretion rates are appropriate for very 
low-luminosity AGN such as Sgr A* in the Galactic 
Center, massive black holes in early type galaxies (e.g., 
iDi Matteo et al.ll2000| ). and some X-ray binaries in qui- 
escence. For higher accretion rates, the electrons cool 
more rapidly than energy can be redistributed by con- 
duction, and so electron conduction is unimportant. Ion 
heat conduction could in principle be relevant at all ac- 
cretion rates since the ion cooling time in RIAFs is al- 
ways longer than the inflow time (by construction); this 
is, however, harder to quantify because the ion conduc- 
tion time is at best comparable to the inflow time of the 
plasma. 

All of the calculations in this paper assume spher- 
ical symmetry. This is formally appropriate only in 
the absence of significant angular momentum, when 
the circularization radius of the inflowing material is 
small compared to the gravitational sphere of influence 
of the central object. We suspect, however, that the 
calculations presented here are relevant to rotating ac- 
cretion flows as well, given the importance of energy 
transport in determining the structure of RIAFs (e.g., 
Blandfo rdfc Begelmanl Il999f) . In a rotating accretion 
disk, however, the conducted energy may be transported 
to the surface of the disk, rather than to large radii. 
This conductive flux might contribute to driving a ther- 
mal outflow. Inde e d, the two-dimensional calculations of 
iTanaka fe Menoul (|2006j ) indicate that conduction helps 
drive a bipolar outflow. Numerical simulations of RIAFs 
with heat conduction are necessary to further quantify 
the importance of conduction and to assess whether any 
resulting heat flux is primarily radial or vertical. Exist- 
ing simulations of RIAFs - which do not include thermal 
conduction - find density profiles with p » 1/2 — 1, in 
which case our calculations indicate that conduction is 
unlikely to lead to super- virial temperatures which would 
be dynamically important (see Figs. [T][3|). Even in this 
case, however, conduction would be important for de- 
termining the electron temperature profile and thus the 
observed radiation. 

Throughout this paper we have treated the dimension- 
less conductivity a c as a free parameter in order to ex- 
plore the effects of conduction in RIAFs. A precise de- 
termination of a c in this context requires understand- 
ing the effective conductivity in a magnetized, turbulent 
plasma. Wave-particle scattering by high frequency tur- 
bulence can limit the mean free path of particles in low 
collisionality plasmas, and thus the conductivity. In the 
Appendix, we argue that instabilities generated by the 
nonlinear evolution of the MRI in RIAFs lead to an up- 
per limit on the ion conductivity of a c < 0.1, in which 
case ion conductivity is unlikely to be dynamically im- 
portant. We find that there is no corresponding limit 
on the electron conductivity, although a large heat flux 
itself may be self-limiting by generating whistler insta- 
bilities (e.g., Pistinner & Eichler 1998). Independent of 
wave-particle scattering, electron transport perpendicu- 
lar to a static magnetic field is suppressed due to the 
small value of the electron Larmor radius. However, 



the rapid separation of neighboring field lines in a tur- 
bulent medium results in an effective isotropic conduc- 
tivity comparable to the field-free conductivity (fjokipiil 
19731 iNaravan fc Medvedevll200H[Cho fc Lazarianl l2004l: 
Chandran fc Maroidl2004t iMaron et alJl2004t ). This sug- 



gests that electron heat conduction will by dynamically 
important even in the presence of a tangled magnetic 
field, with a c ~ 0.1 — 1. With this value of the elec- 
tron conductivity, we find that conduction significantly 
modifies the dyn amics of Bondi accretion (see Figs. [7J[8|). 

iBalbusI (|2000l ) has shown that low collisionality con- 
ducting plasmas are unstable to convective-like motions 
in the presence of an outwardly decreasing tempera- 
ture gradient (the MT I). In a non-rotating atmosphere, 
IParrish fc Stond (|2005l ) demonstrate that the MTI leads 
to magnetic field amplification and a substantial heat flux 
down the temperature gradient. It is natural to specu- 
late that our solutions represent the non-linear saturation 
of the MTI in accretion disks. However, simulations of 
the MTI in RIAFs are necessary to test this hypothe- 
sis since it is not clear what effect di fferential rotation 
and th e MRI will have on the results of IParrish fc Stonel 
(2005) . The MTI will also be important in purely spher- 
ical accretion, as indicated by Figure^ which shows the 
growth rate of the MTI in units of the inverse infall time 
calculated from one of our spherical accretion solutions 
with / e = and a c = 0.1. As Figure [9] shows, the 
MTI will be particularly important at large radii and 
will likely significantly modify existing results on the im- 
p act of magnetic fields on spherica l accretion (e.g., those 
of llgumenshchev fc Naravanll2002l ). 

Although our focus in this paper is on the electron 
energetics in RIAFs, we have also carried out some pre- 
liminary calculations including the ions (neglecting ion 
conduction for reasons discussed above and in For 
spherical accretion with significant ion heating, the ac- 
cretion rate is reduced relative to the Bondi rate by even 
more than is shown in our electron-only calculations in 
Figure [8] This is consistent with the arguments in 
The ions are non-relativistic throughout most of the flow, 
and therefore have 7* ~ 5/3. In the presence of ion heat- 
ing, the effective adiabatic index of the flow is > 5/3, in 
which case there is no transonic Bondi solution in New- 
tonian gravity. Instead, the inflowing matter is in ap- 
proximate hydrostatic equilibrium throughout the bulk 
of the flow and the solution never achieves the free-fall 
density profile of p oc r~ 3 / 2 (see the discussion around 
e q- |21j). The further reduction in accretion rate due to 
the inclusion of ion heating (even absent ion conduction) , 
while strengthening our conclusions, is a consequence of 
the singular nature of the 7^ ~ 5/3 Bondi solution, for 
which Rs — 0. Including additional physical effects such 
as rotation will partially remove this singular behavior 
and may be required to realistically assess the impact 
of ion heating and conduction on the dynamics of hot 
accretion flows. 

As noted in the Introduction, radiation emitted from 
close to the accreting object can play a role analo- 
gous to that of conduction in heating up the outer re- 
gion of the accretion flow to super-virial temperatures. 
This me chanism has been studi e d in the con text of RI- 
AFs by iPark fc Ostrikerl lfl998l H999l [20f)l . The ac- 
cretion flow in that case can exhibit time dependence 
due to the lack of a steady-state solution or to ther- 
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mal instability (lOstriker et all Il976t ICowie et all 119781: 
Bisnovatvi-Kogan fc Blinnikovl Tl980l : iKrolik fc Londonl 



1983). We do not expect time dependence in our prob- 
lem, since our calculations show that a steady-state solu- 
tion does exi st and conduction tends to stabilize thermal 
instabilities (jKrolik fc Londonl 119831 ): a time-dependent 
calculation or a stability analysis would, however, be re- 
quired to confirm this. We also note that our calculations 
are only relevant for systems with sufficiently low lumi- 
nosities as to be unaffected by radiative preheating (as 
discussed above, our solutions require low luminosities so 



that a significant fraction of the accretion energy can be 
conducted outwards rather than being radiated away). 
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APPENDIX 
AN UPPER LIMIT ON a c 

In a collisionless plasma, the temperature and pressure need not be the same perpendicular (Tj_) and parallel (Tn) 
to the local magnetic field. The temperature anisotropy cannot, however, be arbitrarily large because high frequency 
waves and kinetic microinstabilities feed on the free energy in the temperature anisotropy, effectively providing an 
enhanced rate of collisions that isotropizes the pressure tensor via pitch angle scattering. This in turn limits the m ean 
free-path of particles and thus the heat conductivity. These issues are discussed in detail in Sh arma et alJ (|2006l ) in 
the context of nonlinear simulations of the MM in a collisionless plasma. Here wc use this effect to estimate an upper 
limit on a c for ions and electrons in RIAFs. For concreteness we focus on instabilities generated by T± > Tn. For 
conditions appropriate to RIAFs, the most important instability regulating the ion (electron) pressure anisotropy is 
the ion cyclotron (whistler) instability. 4 

For electromagnetic fluctuations slow compared to the cyclotron frequency, the magnetic moment (fi oc T±/B) 
is an adiabatic invariant. As a result, temperature anisotropies with T± =/= Tji are created due to the fluctuating 
magnetic field associated with the MRI. Forma lly, a plasm a with any nonzero temperature anisotropy can be unstable 
to the ion cyclotron and whistler instabilities (jStbd 11992 ). However, there is an effective anisotropy threshold given 
by the requirement that the unstable modes grow on a timescale comparable to the disk rotation period. Gary and 
collaborato r s hav e ana lyzed the ion cyclot ron instability in detail through linear analysis and numerical simulations. 
iGarv et ail (|1997l ) and iGarv fc Led (|1994f ) calculate the anisotropy required for a given growth rate 7 relative to the 
ion cyclotron frequency fij 

--*>#- (AD 

where = P\\,i/(B 2 /8n) and where S — 0.35 and q — 0.42 are fitting parameters quoted in equation (2) of lGarv fe Led 
(1994) for 'j/fli — 10~ 4 . Moreover, for 7 C fii the threshold anisotropy depends only very weakly on the growth 
rate 7. As a result, equation (|A1[) provides a reasonable estimate of the pressure anisotropy required for pitch angle 
scattering by the ion cyclotron instability to be important on a timescale comparable to the disk rotation period (il -1 ). 

If the characteristic timescale for the magnetic field to change by order unity due to turbulence in the accretion 
disk is then the rate of pitch angle scattering (the "collision frequency") required to ensure that the pressure 

anisotropy never exceeds the threshold given by equation (|A1[) is v5p w Hp — > v w £1/3^ J S where Sp — p± — p\\ is the 
pressure anisotropy at the threshold in equation (IA1|) . For pitch angle scattering at a rate u, the ion conductivity is 
given by 

■ 2 'H\ ( _S 

PI 



Ki « _L « nVl r ( - ) ( " (A2) 



where the scale- height of the disk is given by H « Uj/fi. Equation (|A2|) corresponds to a c ~ S/3,, \H/r. Nonlinear 

simulations of the MRI in RIAFs imply f3 ~ 10 in the bulk of the disk, in which case a c w O.lH/r for the ions. This 
is, of course, only an upper limit because there may be additional sources of wave-particle scattering and because 
equation (|A2|) is the conductivity along field lines while the conductivity used throughout the rest of this paper is an 
effective isotropic conductivity. 

Th e argument for an upper limit on the electron conductivity proceeds analogous to that above, with equation 
(jAll) being replaced by t he correspondi n g lim it on the electron anisotropy induced by the whistler instability. For 
non-relativistic electrons. IGarv fc Wang (1996) show that the threshold for the whistler instability is similar to that 



of equation (|A1[) with /3u j — » /3\\. e - Thus the electron conductivity is given by 

Equation (|A3j) is far less restrictive than equation (|A2j) . generally implying a c > 1 for the electrons, i.e., no limit on 
a c at all (saturation of the electron heat flux, however, implies a c < 1). This conclusion strictly only applies when 
the electrons are non-relativistic and more work is needed to understand the extension of this argument to relativistic 
electrons. 



4 The mirror instability provides a comparable rate of ion pitch angle scattering ( Sharma "et~aTl l2006). 



9 




j I i I i I i I i I 

10 10 2 10 3 10" 10 5 



r (Rs) 

Fig. 1. — Electron temperature for a fixed density profile with n oc r~ p , with p varying from 0.5 to 1.5 in increments of 0.1 (from bottom 
to top). The dotted line is the virial temperature T„; r oc r _1 . Values for the other parameters are a c /a = 1 and f e = 0.1. For p > 1, these 
solutions are unphysical because the electron temperature is super-virial. 
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Fig. 2. — Electron temperature with n oc r °' 5 for various a c /ot and / e . The thick solid line is the calculation from Figurc[T]with p = 0.5, 
a c /a = 1 and f e = 0.1. The solid curves correspond to a c /a varying from 0.01 to 1000 by factors of 10 (from top to bottom). The dashed 
curves are for f e = 0.5 (top), 0.01 and (bottom). The dotted line is the virial temperature T v i r oc r —1 . Note that for this density profile 
the solutions are never super-virial. 
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Fig. 3. — Electron temperature with n oc r 1 for various a c /a and / e . The thick solid line is the calculation from Figure HI with p = 1, 
a c /a = 1 and f e = 0.1. The solid curves correspond to a c /a varying from 0.1 to 1000 by factors of 10 (from top to bottom). The dashed 
curves are for f e = 0.5 (top), 0.01 and (bottom). The dotted line is the virial temperature T v i r oc r _1 . 



12 



Johnson & Quataert 




Fig. 4. — Radial temperature profile for the solution of the Bondi accretion problem with heating and conduction, for f e = 0.3. The 
solid curves correspond to a c varying from 1 to 10 4 by factors of 10 (from top to bottom). The dashed curves are for a c = 10~ 2 (top) and 
0.1 (bottom). 
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Fig. 5. — Radial density profile for the solution of the Bondi accretion problem with heating and conduction, for f e = 0.3. The solid 
curves correspond to a c varying from 1 to 10 4 by factors of 10 (from bottom to top). The dashed curves are for a c = 10~ 2 (top) and 0.1 
(bottom). 
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Fig. 6. — Radial velocity profile for the solution of the Bondi accretion problem with heating and conduction, for f e = 0.3. The solid 
curves correspond to a c varying from 1 to 10 4 by factors of 10 (from bottom to top). The dashed curves arc for a c = 10~ 2 (top) and 0.1 
(bottom). 
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Fig. 7. — Sonic radius (in units of the Schwarzschild radius R,s) as a function of a c for f e = (top), 0.03, 0.1, 0.3 and 0.5 (bottom). 
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Fig. 9. — MTI growth rate in units of the inverse infall time f or a spherical accretion solution from [|4] f e = and a c = 0.1. Plotted as a 
function of radius is the maximum growth rate (equation [29] of Balbus 2000), achieved when the conduction time is short but the tension 
forces are small. 



